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We analyze in mathematical detail, within the framework of the QMUPL model of spon- 
taneous wave function collapse, the von Neumann measurement scheme for the measurement 
of a 1/2 spin particle. We prove that, according to the equation of the model: i) throughout 
the whole measurement process, the pointer of the measuring device is always perfectly well 
localized in space; ii) the probabilities for the possible outcomes are distributed in agreement 
with the Born probability rule; iii) at the end of the measurement the state of the micro- 
scopic system has collapsed to the eigenstate corresponding to the measured eigenvalue. This 
analysis shows rigorously how dynamical reduction models provide a consistent solution to 
the measurement problem of quantum mechanics. 



I. INTRODUCTION 

In standard textbooks on Quantum Mechanics, e.g. in [1], one can find the following axioms 
defining the quantum theory: 

Axiom 1: states. A Hilbert TC space is associated to each physical system and the state of the 
system is represented by a vector in 7i. (In the following, we will always assume vectors to be 
normalized.) 



Axiom 2: observables. To any observable quantity of the system is associated a self-adjoint 
operator in 7i. The only possible outcomes of a measurement of an observable are the eingenvalues 
of the associated operator. 
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Axiom 3: Schrodinger equation. Given l-^o) the state of the system at an initial time to = 0, 
its state at any subsequent time t is represented by \tpt), which is the solution of the Schrodinger 
equation: 

zn^lVt) = H\^Pt), (1) 

for the given initial condition; the self-adjoint operator H is the Hamiltonian operator associated 
to the system. 

Axiom 4: Born rule. Let l'^) be the vector describing the state of the system at a given time; 
then the probability that the outcome of a measurement of an observable A at that time is one of 
the values a„ belonging to the spectrum of A, is given by the Born probability rule: 

P[a„] = {il^lPnm, (2) 

where P„ is the projection operator associated to the eigenmanifold of the operator A correspond- 
ing to the eigenvalue a„. 



Axiom 5: wave-packet reduction. At the end of a measurement process the state of the 
system changes according to the rule: 

^..^—.^ mp^' 



where P„ is the projection operator associated to the outcome a„ of the measurement. 

As well known, the last axiom gives rise to the measurement problem in Quantum Mechanics, 
because of which the theory, as it stands, cannot be considered a consistent description of physical 
phenomena. Many tentative solutions have been suggested, among which dynamical reduction 

models are one of the few promising proposals; their general structure has been already fully 
described in the past literature [2, 3]; here we limit ourselves to list the axioms defining them (at 
the non relativistic level): 



Axiom A: states. A Hilbert H space is associated to each physical system and the state of the 
system is represented by a (normalized) vector in H. 
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Axiom B: evolution (continuous version). Given the system initially in a state described by 
the vector |^o)) its state at any subsequent time t is represented by l-^t), which solves the following 
stochastically modified Schrodinger equation: 



-'-Hdt + V\{A-{A)t)dWt - ^{A-{A)tfdt 



(4) 



where Wt is a standard Wiener process defined on a probability space {0,,J^,F), while 
{A)t = {iptlAlijjt) is the quantum average value of the operator A, which is a suitably chosen 
(in this case, self-adjoint) operator; A is a positive constant controlling the strength of the collapse^. 



Axiom C: ontology. Let V'(xi,X2, . . .xjv) = (xi,X2, . . .xjvIV') the wave function for a system 
of N particles (which for simplicity we take to be scalar) in configuration space. Then 

n[''\y:n) = rUn J d^xi . . . . . . d^XAr |V'(xi,X2, . . . xjv)P (5) 

represents the density of mass^ of the n-th. particle of the system, to which a total mass m„ is 
associated [4, 5]. 

Axiom A is equal to 1, while B replaces 3. Indeed, B embodies 3, meaning with this that a 
sensible choice for A and A can be made such that Eq. (4) practically reduces to Eq. (1) when a 
microscopic quantum system is taken into account: sec [6] for a recent and exhaustive review of 

the subject. 

The remarkable property of collapse models is that also the other axioms of Quantum Mechanics 
derive from B (and C, of course): the aim of this paper is to show how axioms 4 and 5 derive from 
axiom B, while in a future paper we will discuss how axiom 2 also derives from B. To be more 
precise, following the previous work of Ref. [7] , we here analyze in mathematical detail, within the 
framework of a specific dynamical reduction model, a von Neumann type of measurement scheme, 
in which a microscopic system interacts with a macroscopic apparatus devised in such a way to 
measure one or more properties of the micro-system. We will show, giving also precise estimates, 
that: 

1. whichever the initial state of the microscopic system, throughout the entire measurement 

^ Eq. (4) can be generalized in different directions [2, 3]; however, its general structure has to be preserved in order 

for the model to provide a solution to the measurement problem. 
^ In the subsequent sections, for simplicity's sake, we will not make reference to the mass density function anymore, 

but we will only keep track of the evolution of the wave function; however it should bo clear that, in order to 

be fully rigorous, all statements about the properties of physical systems should be phrased in terms of their 

mass-density distribution, not in terms of the wave function. 
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process the apparatus has a definite position in space, its wave function being always ex- 
tremely well localized; 

2. the only possible outcomes correspond to those given by standard quantum mechanics, with 
probability almost equal to 1; 

3. the probability of getting a certain outcome is given by the Born probability rule within an 
exceedingly high degree of approximation; 

4. after the measurement, the state vector of the microscopic system collapses to a state which 
practically coincides with the eigenstate of the measured observable, corresponding to the 
eigenvalue which has been observed. 

Needless to say, these properties were already known since very long time, and indeed they represent 
the very motivation behind the original GRW model [8] and its subsequent generalizations, and 
the reason for its success; our goal here is to derive them in a rigorous mathematical way from the 
equations of a specific model of wave function collapse. 

The paper is organized as follows. In section II we will introduce the measurement model our 
analysis is based upon, and we will discuss its physical features. In section III we will study the 
special case in which the microscopic system has been prepared in an eigenstate of the operator 
associated to the observable the model is devised to measure, while in section IV we will analyze 
in full detail the case of an arbitrary initial state. In the concluding section V we will summarize 
the features of our model and draw our final conclusions. 

II. THE MEASUREMENT MODEL 

We begin our discussion by presenting the measurement model we will use in the following 
sections: the setup consists of a microscopic system S interacting with a macroscopic system A 
which acts like a measuring apparatus; both systems are described in quantum mechanical terms. 
Here below we give the details. 

A. The microscopic system 

We consider a single measurement process, in which the experimenter is able to distinguish 
among a finite set of outcomes. Accordingly, we assume that the microscopic system S can be 
described, for what concerns the measurement process, by a finite-dimensional complex Hilbert 
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space. For the sake of simplicity, and without loss of generality, we can consider the simplest case: 
Hs = C^, because the generalization of what follows to is quite straightforward. Since the most 
general self-adjoint operator O acting on can be written as 



where |+) and |— ) are the eigenstates of O, while o+ and o_ are its two real eigenvalues, for 
definiteness and with no loss of generality, in what follows we will take o± = ±h/2 and O to be 
the z-component of the spin, Sz, of a 1/2 spin particle. 



We take the following model for the measuring apparatus A, which is general enough to describe 
all interesting physical situations: we assume that the apparatus consists of a fixed part plus a 
pointer moving along a graduate scale, in such a way that different positions of the pointer along the 
scale correspond to different possible outcomes of the measurement. To simplify the analysis, we 
study the evolution of the center of mass of the pointer only, and disregard all other macroscopic 
and microscopic degrees of freedom; accordingly, the pointer will be treated like a macroscopic 
quantum particle of mass m moving in one dimension only, whose state space is described by the 
Hilbert space Ha = L2(M). 



We assume that the pointer of A undergoes a spontaneous collapse mechanism according to the 
Quantum Mechanics with Universal Position Localization (QMUPL) model first introduced in [9] 
and subsequently analyzed in [10] (see also references therein), while the microscopic system S 
evolves according to the standard Schrodinger equation, since, as typical of dynamical reduction 
models, the stochastic collapse terms have very little effects on microscopic quantum systems. 
Accordingly, we take for the evolution equation of the composite S + A system the following 
stochastic differential equation^ (SDE) defined in the Hilbert space H = Hs® Ha '■= (8) L^(R): 



= o+|+)(+| + o_|-)(- 



(6) 



B. The measuring apparatus 



C. The Dynamics 



d\^t)= -'^Htdt + VXiq-{q)t)dWt-^iq-{q)tfdt , 



(7) 



^ See [11] and [12] for theorems on the existence and uniqueness of solution for this type of equation. 
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which is precisely of the form (4) , where Hf is the (time dependent) standard Hamiltonian operator 
of the composite system, and q is the position operator associated to the centre of mass of the 
pointer.^ In the following we will use capital letters (^', . . .) to denote a state vector for the 
composite S + A system, and lower case letters {tp, (j),. . .) to denote a state vector referring to the 
pointer alone. 

As discussed in [10] we take for the constant A appearing in (7) 

777 

A:^— Ao, (8) 
mo 

with mo — 1.7 x 10~^^ Kg being a reference mass of order of a nucleon mass and Aq ^ 10~^ m"^ 
sec~^. For definiteness, let us consider a pointer of mass m = 1 g (i.e., a pointer made of an 
Avogadro number of nucleons), and let us define, for later convenience, the quantities 



hx h 

Lo:=2\ — ~ 5.0 X 10-^ sec-i and (t„ := W ~ 4.6xl0-^^m. (9) 

V m V mu 

We take the Hamiltonian Ht to be of the form Ht = Hs+H^+Hj^t- The first term is the quan- 
tum Hamiltonian for the microscopic system: we assume that the time scale of the free evolution of 
the microscopic system is much larger than the characteristic time scale of the experiment ( "instan- 
taneous measurement" assumption); accordingly we take Hs to be the null operator. The second 

term is the quantum Hamiltonian of the pointer, which wc take equal to that of a non-relativistic 
free quantum particle of mass m: H^^ = p'^ /{2m), where p is the momentum operator. 

Finally, we assume the interaction term i?iNT between the two systems to be of the von Neumann 
type, and devised in such a way to measure the spin operator Sz'- 

HiMt) = KAjSz^p, (10) 

where k is a coupling constant and A-^: t i— > is a T-normalized,^ non negative, real valued, 
function of time, identically equal to zero outside a given interval of the form (to, to + T), i.e., 
outside the time interval of length T, say T = 1 sec, during which the experiment takes place; we 
choose the time origin in such a way that the experiment begins at io = sec. As it is well known. 



* Thus, strictly speaking, we should write Is ^ q for the position operator for the pointer, where Is is the identity 
operator in Hs ■ We avoid such a way of writing, when no confusion arises. 



By a T-normalized function, we just mean 



Afdt= / Ajdt = T. 

J to 



Note that Af depends also on the initial time to; we will avoid to indicate it explicitly, when non confusion arises. 
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-ffiNT generates the following type of evolution, depending on the initial state of the micro-system 
S: 

[c+ 1+) + c_ I-)] ® |<^o) ^ c+ 1+) ® \(f)+) + c_ I-) ® |<^_) , (11) 

where \(p±) are final pointer states spatially translated with respect to the initial state |^o) by the 
quantity ±{h/2)KT. 

The strength of the coupling constant k has to be chosen in such a way that the distance hnT 
between the initial state |<^o) of the pointer and any of the two final states is macroscopic; for 
definiteness, let us choose hK = l cm sec~^, so that HkT = 1 cm. 



D. The initial state 



We take the initial states of the microscopic system S and of the macroscopic apparatus A to be 
completely uncorrelated, as it is customary and appropriate for the description of a measurement 
process. Accordingly, we assume the initial state of the total system 5 + .4 to be 

[c+|+) + c_|-)]®|</.o), (12) 

where |0o) describes the "ready" state of the macroscopic apparatus A. 

Regarding the initial state \(I)q-) of the pointer, some considerations have to be done. In [10] 
it has been shown that, according to equation (7), the wave function for the centre of mass of an 
isolated quantum system reaches asymptotically (and very rapidly, for a macro-object) a Gaussian 
state of the form 




C(^) = WtT— 2 exp 



-^^{x - xtY + iktx 



(13) 



(modulo a time-dependent global phase factor) with Gq defined as in Eq. (9). For later reference, 
let us observe that the dispersion of the Gaussian function of Eq. (13) in momentum space is 

Gp = ~ 1.6 X 10~^^ kg m sec-\ (14) 

quite close to the minimum allowed by Heisenberg's uncertainty relation, and that the centres of (f)f 
in position and momentum space are given by xt and hkt respectively. 

In our measurement model, we assume that the pointer is isolated for the time prior to the 
experiment; during this time, as shown in the past literature, its wave function converges rapidly 



8 



towards a state close to (13), which we therefore assume to be the initial state of the pointer. To 
summarize, we take as the initial state of the composite system S + A the ket 

|^'o) = [c+|+) + c_|-)]®|G,0), (15) 

where (x|G,0) is of the form (13). We choose the natural reference frame where the pointer is 
initially at rest, so that kg = m~^, with the origin set up in such a way that xq = m. 

III. MEASUREMENT OF AN EIGENSTATE 

We begin our study of the model by looking for the solution of equation (7) satisfying the initial 
condition 

|*±) = |±)0|G,O), (16) 

where the symbol it means that the state |ib) is either |+) or |— ), i.e., an eigenstate of the 
operator Sz- We will show that, in this special case, the state of the microscopic system does 
not change in time, while the pointer moves along the scale so to give the correct outcome of the 
measurement. 

A. The linear equation 

Following the standard procedure outlined e.g. in [10] we pass from the non-linear equation (7) 
to the corresponding linear equation 

^t) , (17) 

where is a standard Wiener process defined on the probability space {0,,J^,Q), the measure 
Q being a new probability measure chosen in such a way that the old probability measure P is 
generated from Q by the martingale The Wiener process of equation (17) is related to 

the Wiener process Wt of equation (7) via Girsonov's rule [13]: 

dWt = d^t-'iV\{q)tdt. (18) 

It is easy to prove that a vector of the form 

|$±) = |±)®|0±> (19) 



Htdt + ^/Xqd^t - t; q'^ dt 
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cpfix). (20) 



solves equation (17), for the initial condition (16), if the wave function (f>f{x) := {x\(f)f) solves the 
following linear SDE, which involves the apparatus degrees of freedom alone: 

B. The solution and its properties 

The solution of equation (20) for the given initial condition is the Gaussian wave function 

(t>f{x) = exp [-at{x - xff + i~kf x + + i Of] , (21) 

whose parameters at G C, and xf ,kf , G M (of obvious meaning) satisfy the following system 
of SDE [10]: 

dat = {X-—aj \ dt (22) 
dxf = [^~kt±^n Af ) dt + ^{d^t- 2VXxt dt} (23) 
dkf = -^/X^\d(t-2^/Xxtdt^ (24) 
djt = (a(x±)2 + lal + ^^dt + yfXxf [d^t - 2VA xf dt] (25) 

= (-2^ - ^ + 4^ ^ T ~^^) + ^1 {^^* - ^4 ' 

where we have denoted by z^' (z^) the real (imaginary) part of the complex number z. 

Eqs. (25) and (26) are of no particular interest in this simple situation, because they just describe 
the time evolution of the irrelevant norm and global phase of the Gaussian solution. Eq. (22) is 
independent of Eqs. (23)-(26); it is deterministic, and easily solved by separation of variables: 

at = tanh ( + cq] , (27) 



4^1 

where Cq sets the initial condition. Eq. (27) determines the time evolution of the spread in position 
and momentum of the Gaussian wave function. In our case, given the initial condition = 
(1 — i)/4(Tg, it follows that at = (1 — i)/4o"^ for all times (i.e., we can set af = l/4cr^ and 
al = —1/4(7^ in equations (23)-(26)): as expected, the position and momentum spreads of the 
wave function of the pointer do not change in time, and remain identically equal to, respectively, 
aq and Up. 

For what concerns equations (23) and (24), which do not depend on equations (25) and (26), 
their solution describes the time evolution of the mean value in position and momentum of the 
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Gaussian wave function. We have to characterize the stochastic properties of the solution with 
respect to the physical probability measure P, i.e. we need to go back to the original noise Wt 
via Girsanov's rule (18), which in this case is very easy, because for a wave function like (21) we 
simply have {q)^ = xf for all times, so that all we have to do is to write dWt in place of the curly 
braces {. . . } in equations (23) and (24): 

dxf = kf ± ^ Af ^ dt + aq yfZ)dWt (28) 
dkf = ^^dWt, (29) 

where we have also taken into account that at = {1 — i)/4:ag for any t > 0. Let us call |^^) the 
normalized physical solutions: l^*^) = l^^t) /II 1*^*^) II' with |<&^) given by Eq. (19); taking also into 
account that xq = m, and ko = m~^, we find the following results. 

1. According to Eq. (29), the average value of the peak of the Gaussian wave function in 

momentum space, (p)^ := (\l'^|p|\I'^) = hk^, does not evolve in time: 

^A{p)f] = = Kg m sec-i. (30) 

2. By Eqs. (28) and (29), the average value of the peak in position space {q)f := = 
xf of the Gaussian wave function evolves in time according to 

= ±^j\j,dt'. (31) 

Equation (31) shows that, when the measurement begins, Ep[(gf^)] moves towards the right 
or the left according to the initial state of the microscopic system S, reaching the final value 
(at the end of the measurement) 

Ep[(g)^±] = = ±0.5cm fort>T. (32) 

3. The variance Vp[(g)^] = Ep[(g)^ — Ep[(g')^]]^ associated to the motion of {q)f is equal to 
the variance computed in [10] (see Sec. VII B), which, for ctj = (1 — i)/4a'^, is given by^: 



±1 ^2 
t \ = 



(ujty iiot) 

cot + ^-^ + ^ ^ 



9' 

2 



(33) 



2 12 

with our choices for the parameters, we have Vp[((7)^] < Vp[(g)^] ~ 1.1 x 10~^^ m^, for any 
t < T. 



We correct in this way a typo contained in Eq. (93) of ref. [10]. 
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Prom the above results we can derive the following important conclusions: 

• Due to the smallness of its variance, the motion of the peak {q)f of the Gaussian wave 
function for the cm. of the pointer is practically deterministic and equivalent to the motion 
of Ep[(g)^], the fluctuations around the mean being so tiny that they can be safely ignored. 
E.g., the probability for {q)f to lie outside an interval of width A centered in Ep[(g)^] can 
be estimated by using Cebicev's inequality; for A = 10~^ cm, we have 

P[| - Ep[(g)±]| > A/2] < 4 YlM)l ^ 4.2 X 10-1^ for any t < T, (34) 

a vanishingly small probability. 

• As such, and because of Eq. (31), the peak {q)f evolves in time as follows: 



± — - Affdt' t<T (+ negligible fluctuations) 

2 Jo 



± — - — = ±0.5 cm t>T (+ negligible fluctuations) 



(35) 



This means that, according to the initial state of the micro-system, the pointer moves in a 
practically deterministic way either towards the left or towards the right, with respect to the 
initial ready-state, displaying in this way the outcome of the measurement. 

• During the measurement, the state of the micro-system does not change. 

This is precisely the expected behavior both for the microscopic system as well as for the macro- 
scopic pointer, when the initial state is given by (19), for an ideal measurement scheme as the one 
of von Neumann here analyzed. 

IV. MEASUREMENT OF A SUPERPOSITION 

Let us now consider the general case where the initial state \si) of the microscopic system S is 
not an eigenstate of Sz, but a superposition of eigenstates of the form 

|si)=c+|+) + c_|-) (|c+|2 + |c_p = l); (36) 

the global initial condition for the micro-system and the apparatus then is 

|*o) = [c+|+) + c_|-)]®|G,0). (37) 
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As in the preceding subsection, we first solve the linear equation, and next move to the non-linear 
one. Due to the linearity of equation (17), its solution, with the given initial condition (37), is 

l^t) = 1+)® + (38) 

where the wave functions and |</>7), in the position representation, are of the form (21) and 
the parameters at,xf solve Eqs. (22) to (26), with the obvious choice of sign and with 

initial conditions: 

\ — i 

"0 = -TT^ = m, A;^ = m'^ 7^ = In |c±|, 6^ = Arg[c±] (39) 

(of course we now assume that c± 7^ 0). 

Since the time evolution of the parameters at is governed by equation (22) which, as we have 
already remarked, is deterministic and decoupled from the evolution equations (23)-(26) for the 
remaining parameters, we observe first of all that the evolution of the spreads in position and 

momentum of the two Gaussian functions \(t)f) and does not change with respect to the case 
analyzed in the previous section. Accordingly, we have = (1 — i)/4c7g for all times, so that the 
spreads of the two wave functions do not evolve and remain identically equal to the asymptotic 
values (jq and a p. 

A. The deterministic evolution of the distances in position and momentum between the 

two Gaussicin components 

Contrary to the preceding case, moving from the solution of the linear equation to the solution of 
the non-linear one is not immediate, since Girsanov's rule (18) involves the quantum average {q)^, 
which in this case turns out not to be a trivial function of the parameters controlling the two 

Gaussian components; namely, one finds that^: 

e 't + e '* 

Of course, this is an entirely expected difficulty, due to the essential non-linearity inherent to 
collapse models; to proceed in the analysis of the problem, it is convenient first of all to analyze 
the evolution of the distance between the maxima of the two Gaussian functions \(f)f) and \4>^), 
both in position as well as in momentum space, and subsequently of their relative weights. 



In Eq. (40) no contribution comes from the overlapping between the two Gaussian components, since each com- 
ponent is coupled to one of the two orthogonal spin state |±), which make the "off-diagonal" terms of the scalar 
product vanish. 
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Following the path outlined in [10], let us consider the differences Xt := — and Kt := 
k^' — k^ , which express at each instant the distance in position and (modulus h) momentum space 
between the centres of the two Gaussian functions and I^om equations (23) and (24), 



keeping in mind that in our case we have af = l/4cr^ and 
following deterministic system for Xt and Kf. 



dt 



Xt 




—u h/m 




Xt 














+ 




Kt 




-2A 




Kt 








-l/4(Tg for all times, we get the 



(41) 



since it does not depend on the noise, it is insensitive to the change of measure and holds true also 
for the non- linear Eq. (7). 

The solution of the above system depends of course on the specific choice for the function A^; 
a simple reasonable choice is the following: 

' 1 t e [0,T] 

(42) 

else, 

which, according to Eq. (35), means that, with restriction to the situation analyzed in the previous 
section, during the measurement the pointer moves at a constant speed either towards the left or 
towards the right, depending on the initial state of the micro-system. According to this choice, Xt, 
given the initial condition Xq = m, evolves in time as follows: 



Xt = < 



-*/^sini^t 



2hK 



-ujt/2 



sin^t-e-^/2sin^(t-T) 



for < t < T, 



for t > T. 



(43) 



Xt ~ 



2 2 

Since lu^^ ~ 2.0 x 10^ sec is a very long time compared to the measurement-time, we can mean- 
ingfully expand Eq. (43) to first order in ut: 

hKt for < t < T = {Hk)-'^ = 1 sec, 

- - \ ) (44) 

1 cm for T < t < uj-'^ ~ 2.0 x 10^ sec 

As we see, the distance between the two peaks increases almost linearly in time, reaching its 
maximum (1 cm) at the end of the measurement process, as predicted by the standard Schrodinger 
equation; after this time, their distance remains practically unaltered for extremely long times, 
and only for t ~ 2.0 x 10^ sec it starts slowly decreasing, eventually going to 0. Note that such a 
behavior, being determined by a;, does not depend on the mass of the pointer, thus a larger pointer 
will not change the situation. The moral is that Xt behaves as if the reduction mechanism were 
not present (as if Aq = 0) so we have to look for the collapse somewhere else. 
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As we shall discuss in the next subsection, the collapse occurs because, in a very short time, 
the weight of one of the two Gaussian wave functions {\(l)f) or \(f)^)) becomes much smaller than 
the weight of the other component; this implies that, when the normalization of the whole state 
is taken into account, one of the two components practically disappears, and only the other one 

survives, the one which sets the outcome of the experiment. Of course, this process is random and, 
as we shall prove, it occurs with a probability almost equivalent to the Born probability rule. 

B. The evolution equation governing the relative weight of the two Gaussian components 

The relative damping between the two Gaussian components of Eq. (38) is measured by the 
stochastic process = 7^ — 7^": if, at a certain time t, it occurs that S> 1, it means that at 
the end of the experiment |^^) is suppressed with respect to 10^"), so that the initial state (37) 
practically evolves to |+) (g) {tp^') (remember that \tpf) = \(j)f) /|| \(f)f) \\ are the normalized states); 
the opposite happens if Ft <C — 1. To be quantitative, let us introduce a conveniently large collapse 
parameter, say 35, and the following definition*^: 

Definition. The superposition (38) is suppressed when |Ft| > 35, i.e., when either 
II 1"?^*^) 11/11 II '^^ reciprocal is greater than e^^ 2± 1.6 x 10-^^. 

Using equation (25) and the Girsonov's transformation (18), we can write the following SDE for F^ 
in terms of the noise Wt associated to the non-linear Eq. (7): 



with initial condition Fq = ln|c-|-/c_|. By using the expression (40) for {q)^, we can re- write the 
above equation as follows: 



This is the result we wanted to arrive at, and we will devote the rest of the section at analyzing its 
physical content. To proceed further with the analysis, it is convenient to perform the following 
time change [14], 



dVt = XXt (2 {q) 



t-xf - x^) dt + \f\Xt dWt 



(45) 



dTt = AX|tanhFtdt + y^XtdWt. 



(46) 



t 




(47) 



The choice made here for the collapse parameter is different from the one made in [10]. We find this new choice, 
which at any rate is arbitrary, more convenient for the problem under study. 
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which allows us to describe the collapse process in terms of the dimensionless quantity s that 
measures its effectiveness. Using Eq. (43), one can solve exactly the above integral and compute 
s as a function of t. Such a function however cannot be inverted in order to get t from s. To this 
end, we use the simplified expression (44) in place of the the exact formula Eq. (43) to compute the 
integral, an expression which, as we have seen, represents a very good approximation to the time 
evolution of Xt throughout the whole time during which the experiment takes place (alternatively, 
we may initially choose in such a way that Xt evolves exactly like in (44) , at least from t = 
to t = T). Accordingly, we have: 

\fc2 2 

s = St 2.0 X 10^^ {t/secf <t <T = 1 sec, (48) 



t = ts ^ \ s ~ (1.7 X 10-^ ^) sec < s < Xh^K^/?, = 2.0 x 10^'^. (49) 

V Xn''-K.^ 

Note that, according to the above equations, the physical time t depends on s through the inverse 
cubic-root of A, i.e. on the inverse cubic-root of the mass of the pointer; this time dependence of 
i on A is important since, as we shall see, it will affect the collapse time. We do not study the 
functional dependence between s and i for t > T since, as we shall soon see and as we expect it to 
be, the collapse occurs at times much smaller than T. 

Written in terms of the new variable s, Eq. (46) reduces to: 

dVs = tanhr^ds + dWs] (50) 

this equation has been already analyzed in [10] , using the theorems of [14] ; here we report the main 
properties. 



1. The Collapse Time 

According to the definition given before, a collapse occurs when \Tt\ > 35; it would seem then 
natural to define the collapse time as the time when \Tt\ first reaches the value 35. However, one 
has to face the event that \Tt\, after having reached such a value, immediately starts decreasing in 

a significant way, jeopardizing in this way the effect of the collapse. To avoid such a possibility, we 
proceed as follows: we will compute the time it takes for |ri| to rcacher a value larger than 35, let 
us say 50, and subsequently we will show that, after having reached such a value, the probability 
that it gets back to a value below 35 is negligible. In this way we can be (almost) sure^ that, once 

^ Here, as well as in the rest of the paper, we use "almost sure" in the physical sense of "with very high probability" , 
not in the mathematical sense of "with the possible exception of a subset of measure 0" . 
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the collapse has occurred, the superposition never re-appears. 

Let us consider the time S = S{uj) when IF^I first reaches the value 50: 

S = inf{s : IF^I > 50}; (51) 

of course we assume that the initial state (36) is such that jFol < 35, otherwise according to our 
definition (as well as for all practical purposes) it would already be a reduced state, not a physically 
meaningful superposition. It can be proven [10] that S is finite with probability 1, and that its 
average value and variance are given by the following expressions: 

Ep[S] = 50tanh50 - rotanhTo, (52) 
Vp[5] = F(50)-F(ro), F{x) = x'^ta-nh^x + xta-nhx-x'^. (53) 

Now, tanh50 ~ 1 — 7.4 x 10"^'^ which is practically 1; let us also consider e.g. the worst case, as 
far as the collapse mechanism is concerned, i.e. the case in which Fq = 0, which means that the 
micro-state is initially in a equal weighted superposition of the two eigenstates. We then have that 
Ep[5] ~ 50 and Y^[S] ~ 50. 

S{u!) is a random variable, so we can not tell exactly when the collapse occurs; since however 
we want to be quite safe that it actually occurs, let us compute the probability that S happens to 
be much greater than, e.g., 10^ times its standard deviation. By a trivial application of Cebicev's 
inequality we have: 



5-Ep[5]| > lOVWl-^] - 7.1 X 10^ 



< 10-^°. (54) 



We can then conclude that, at time s = Ep[S] -|- 10^\/Vp[5] ~ 7.1 x 10^, the collapse has almost 
certainly occurred (with probability greater than 1 — 10~^°) and that it is an irreversible process 
(as we shall soon prove). Moving back from the effective time s to the physical times t by using 
Eq. (48), we then define the collapse time as follows: 



y All K 

the collapse occurs within a time interval smaller than the perception time of a human observer. 
The above formula shows also that, as expected, Tc is proportional to the inverse cubic-root of 
the mass of the pointer (since A is proportional to the mass): the bigger the pointer, the shorter 
the collapse time. With our choice for Aq, even for a 1-g pointer the reduction occurs practically 
instantaneously. 
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It is important to note that, at time Tc ~ 1.5 x 10~^ sec, the distance between the two Gaussian 
components is approximately Xtq — 1.5 x 10""^ cm: this means that, with very high probabihty, 
the collapse occurs before the two components have enough time to spread out in space to form a 
macroscopic superposition. This means that, from the physical point of view, there is no collapse 
of the wave function at all, since it always remains perfectly localized in space at any stage of the 
experiment. In any case, we will keep talking of collapse of the wave function, meaning with it 
simply the event \Tt\ > 35. 



2. The Collapse Probability 

Let us call P+ the probability that Fg hits the point +50 before the point —50, i.e. the probability 
that \(l)f) survives during the collapse process so that the outcome of the measurement is "+^/2". 
Such a probability turns out to be equal to [10]: 

1 tanh 50 + tanh Fq 

while the probability P that Tg hits the point —50 before the point +50, i.e. that the outcome of 
the experiment is "— /i/2", is of course: 

1 tanh 50 — tanh Fq 
" ~ 2 tanh 50 ' ^ ' 

By taking into account that tanh 50 ~ 1 — 7.4 x 10~^^ 2± 1, we can write, with extremely good 
approximation: 

1 e'^o e'^'^o 
P+ ~ - [1 + tanhFo] = -p jT- = — T = |c+p, (58) 

P_ ~ - [1-tanhFo] = ^ ^ = — ^ = |c_P. (59) 

We see that the probability of getting one of the two possible outcomes is practically equivalent 
to the Bom probability rule! On the one hand, this is an entirely expected results, since collapse 
models have been designed precisely in order to solve the measurement problem and in particular to 
reproduce quantum probabilities; on the other hand, it is striking that a very general equation like 
Eq. (7), which is meant to describe both quantum systems as well as macroscopic classical objects 
(i.e. all physical situations, at the non relativistic level), when applied to a measurement situation, 
provides not only a consistent description of the measurement process, but also reproduces quantum 
probabilities with such a good precision. 
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3. Stability of the Collapse Process 

We have already anticipated that, since Vg evolves randomly, there is the chance that, after 
having reached e.g. the value +50, i.e. after that the wave function collapsed to the state |^"^), 
it becomes smaller than 50 instead of keeping increasing, eventually getting closer and closer to 
0. When such an event occurs, the superposition of the two Gaussian wave functions, which 
was previously reduced, reappears again, jeopardizing in this way the entire collapse process and 
localization properties of the pointer. We now give an estimate of the probability for such an event 
to occur. 

Let us call (5+ the probability that F^, after having reached the value +50 at time 5, does not 
go back to a value smaller 35: 

(60) 

such a probability turns out to be [10]: 

Q+ > (l+tanh50)-^^^^^i^4^ ~ 1-9.3 X 10-l^ (61) 
1 + tanh 15 

which is practically equal to 1: once a localization occurs, the superposition can de facto never 
re-appear. 



inf Tg 

s>S 



> 35 



C. State vector after the collapse 

At time t>Tc the normalized sate vector = with given in (38), is: 

1+) \G+,I) + (,]-) \G-.t) 



(62) 



where et = e ^'^^ ^ and the normalized Gaussian states {x\G±,t) are defined as follows: 

1 - i 



{x\G±,t) 



27ra2 



exp 



4a| 



■(x - xff + ikf x + i of 



(63) 



Let us assume that the collapse occurred in favor of the "+^/2" eigenvalue, i.e. in such a way that 
Vt > 35 for t>Tc, with very high probability; it follows that: 



et < e-^^ ~ 6.3 x 10"^^ ^ t>Tc, 



(64) 



and we can write, with excellent accuracy: 



(65) 



19 



We recover in this way the postulate of wave packet reduction of standard quantum mechanics: at 
the end of the measurement process, the state of the micro-system is reduced to the eigenstate 
corresponding to the eigenvalue which has been obtained as the outcome of the measurement, 
the outcome being defined by the surviving Gaussian component {\G+,t) in this case). Note the 
important fact that, according to our model, the collapse acts directly only on the pointer of the 
measuring apparatus, not on the micro-system; however, the combined effect of the collapse plus 
the von Neumann type of interaction is that the microscopic superposition of the spin states of the 
micro-system gets rapidly reduced right after the measurement. 

Note finally that, after the collapse, the states of the micro-system and of the pointer are de facto 
factorized: as such, after the measurement process one can, for all practical purposes, disregard 
the pointer and focus only on the micro-system for future experiments or interactions with other 
systems, as it is custom in laboratories. 



D. The end of the experiment 

In this final subsection we study how, after the collapse, the "winning" component {\G+,t) or 
\G—,t)) moves in space, i.e. how their centers xf or evolve in time, whether they move in such 
a way to display the correct outcome of the measurement. To this purpose let us define: 

Xt = xi + x^, Kt = kt + K, (66) 

so that xf and x^ as functions of Xt and Xt are given by: xf = {Xt+Xt)/'^' and x^ = —{Xt—Xt)/2. 
Prom Eqs. (23) and (24), taking also into account (40), one finds out that Xt and Kt satisfy the 
following SDEs: 

dXt = —ktdt + ojXt tanh dt + 2yJuOg dWt, = m, 

m (67) 

dKt = 2XXttanhrtdt + 2VXdWt, Kq = m-\ 

where Xt is given by Eq. (43). This is a non- linear system, since it depends in a non-linear way on 
Tt, which is also a stochastic process; as such, (to our knowledge) the system can not be exactly 
solved. To circumvent this problem, let us consider the following two auxiliary linear systems: 

h 

dXf" = — dt ± ujXt dt + 2Vwo-„ dWt, = m, 

* m * V g t, > ^gg^ 

dkf = ±2\Xtdt + 2^/\dWt, k^=0 m'^ 

(with an obvious meaning of the symbols), which have been obtained in the first case (-I-) by 
replacing tanhT^ with -|-1, and in the second case (— ) by replacing tanhT^ with —1. Clearly, we 
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have: Xf < Xt < X^ and Kf < Kt < for any t such that > 0, which is true for all the 
time during which the experiment takes place, and much longer. Such linear systems can be easily 
solved; concerning X^, and after some tedious calculations one finds the following time dependence 
for the mean: 



and for the variance: 



-Xt + hKt for t < T, 

(69) 

-Xt + hKT for t > T, 



Vp[X±] = Aal 



12 



(70) 



We use the above results to approximate the time evolution of Xt and thus of x'l and x^, which 
we are interested in. We consider separately the case t <Tc (before the collapse) and t>Tc after 
the collapse: in the first case, we cannot control the behavior of Tt, thus the most we can say is 
that I tanhrd < 1, which has already been used to bound Xt between Xf and X^; in the second 
case, we know that with very high probability | tanhTfl > tanh35, which is a very strong bound. 
Case 1, before the collapse: t < Tc- Within this time interval, the two Gaussian components 
|G+,t) and |G— ,t) start separating, as Xt increases in time; in particular, at time t = Tq, when 
the collapse has (almost certainly) occurred, we have: 

Et?[X^J ~ ± ^ hncvT^ ~ ±5.9 x lO'^^m, (71) 

which has being obtained from Eq. (69) by expanding Xt, as given by Eq. (43), to second order in 
ut; moreover, we have from Eq. (70): 

Yw[X^J ^ Aua^Tc ~ 6.5 x lO'^V^ (72) 

This means that, on a macroscopic scale, X^^ 2± Ep[X^^]; since Xtq — HkTc ^ 1.5 x 10~^ 
m S> Ep[X^^], and keeping in mind that Xi^^ < Xj,^ < X^^, we can write, with very high 
probability and very good approximation: 

~ +^Xtc ^ +\hKTc ^ +7.7xl0-^m, (73) 
~ ~Xtc ^ —huTc c± -7.7xl0'~^m. (74) 

Accordingly, and as expected, the two components move symmetrically in opposite directions, one 
towards the right and the other towards the left, but not fast enough for a macroscopic superposition 
to occur, before the collapse enters into play and suppresses one of them. 
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Case 2, after the collapse: t > Tc- Let us assume that the cohapse is such that the outcome 
"+^/2" is given; this means that almost certainly > 35, V i > Tc- Given this, let us first of all 
show that Xf remains very close to X^', for very long times; then, by approximating X^ with X^, 
we will show how xf and evolve in time. 

From Eqs. (67) and (68), taking into account that tanhF^ > —1, we find: 

rt rt 



k+-kf = 2X f Xt>{l-tanhTt')dt' < AX f Xt> dt' = 
Jo Jo 

~ 2XhKt'^ (75) 



'0 Jo 

8XhK r , rf/9 / i^t . iot 
1 _ e-^*/2 ( cos _ + sin _ 



a;2 



and 



X+-X^ = — [Kj - K^, )dt' + UJ / X*/ (1 - tanh Tf) dt' 
Jo Jo 

[ {k+ - k^, ) dt' + 2u I Xf dt' < -2{Xt - hnt) ~ hnut^ . (76) 
Jo Jo 



< 

m 



At time t = Tc, we then have: - kj.^ ~ 2XhKT^ ~ 2.8 x 10^^ m'^ and: X^t^ - Xj.^ ~ 



HkloT^ ~ 1.2 X 10-^^ m. 



We use these results as initial conditions, at time Tc, to find, by integrating once more Eqs. (67) 
and (68), and by using the two inequalities tanhFj > r) = tanh 35, ^ t > Tc and < ^ 2± 1 cm, 
the following estimates: 

ft 

'Tc 

and 

ft ft 



k+ -kf = k+^ - kxc + 2x1 Xt,{l - tanhF^/) dt' < k+^ - kxc + 2Xrie{t - Tc) (77) 

J Tn 



X^ -Xi- < X+, -Xj,^ + — j {k+ - k^, )dt' + UJ j Xt' (1 - tanh Tf) dt' 

ITT' Jtc JTc 

< X+^ - X^^ + -(K+^ - k^^){t - Tc) + ^rjiit - Tc)^ + coriiit - Tc) 



m 

2 2 

~ loHkT^ + y hKT^{t -Tc) + ^ r]£{t - Tc)^ + ur]l{t - Tc) 

~ (l.2x 10-^^ + 2.9x 10-^'^t + 5.0x 10-^2t2)m, (78) 

We see that for very long times, by far much longer than the time during which the experiment 
takes place, the distance between Xj,^ and X^^ remains small, so small that we can replace 
with Xj,^ for all practical purposes. 

On the other hand, Xj,^ is, on a macroscopic scale, very close to its average value EpfX^"*"], its 
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variance, as given by Eq. (70), being extremely small; accordingly we have: 



Xt + X, 



t Xt + xt 



Xt + MXt 



HkT 



+ 



t < T, 



t > T, 



(79) 



which is the desired result: the pointer, represented in this case by \G+,t), moves at a constant 
speed towards the right and stops at the position hKT/2, displaying in this way the correct outcome. 

To conclude the analysis, let us see what happens also to the other component, \G—,t), which 
has been suppressed by the spontaneous reduction process. Its center moves approximately as 
follows: 



X,-X, 



Xt-xt 



+- 



hnt 
HkT 



t < r, 



(80) 



i.e. the negligible Gaussian component moves to the left of the graduate scale, but then slowly 
converges towards the other wave function. 

As a final remark, we note that, at very long times of order uj~^ ~ 2.0 x 10^ sec, the statistical 
fluctuations become relevant also on the macroscopic scale, thus approximating any actual value 
with its statistical average becomes less and less precise. However, times of order lo~^ ~ 2.0 x 10^ 
sec arc by far much longer that the time required for the experiment to end; moreover, for such 
long times the assumption that the global system is isolated certainly looses its validity; the 
measurement model should then be refined, in order to include so long time scales. 



V. CONCLUSIONS 



In the present work we have analyzed the quantum theory of measurement within the framework 
of dynamical reduction models, resorting to the von Neumann type scheme of measurement process 
and to the QMUPL model of spontaneous wave function collapse. We have proven the properties 
listed in the introductory section, showing in this way how the axioms 4 and 5 of standard Quantum 
Mechanics arise in quite a straightforward way from the dynamical evolution law governing models 
of spontaneous wave function collapse. 
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We hope that our analysis makes clearer the mechanism with which dynamical reduction models 
provide, at least at the non relativistic level, such an accurate description of measurement processes, 
and more generally of all physical situations. 
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APPENDIX A: A MISTAKE CORRECTED IN REF. [10] 



When K = 0, i.e. for a free particle of mass m moving according to the SDE: 

^2 A 



diptix) 



l^dt + VXiq- {q)t) dWt -^{q- {q)tfdt 



the two equations (25) and (26) for and Ot, respectively, become (we neglect the ±): 



A 



dOt = {-^^--af + 
' 2m * m * 



dt + ^/Xxt jd^t - 2VXxt dij 

dt + VX^ xt {c?6 - 2VXxt dt^ , 



4KF 



(Al) 



(A2) 
(A3) 



which differ from the corresponding Eqs. (12) and (13) of [10], in the first case for the extra factor 
A/4af and in the second case for the factor AaJ/4(af )^. We correct in this way a mistake made 
in Ref. [10], which however does not affect the other results contained in that paper. 
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